use "$Data_generated/Diocese", clear
cibar secburial_total if year==1935, over(eveque_antimo) barc(gs12 gs3) gr(ytitle("Average number of Secural Funerals") legend(order(1 "Non Anti-AF Bishops" 2 "Anti-AF Bishops") pos(6)))
graph export "$Figures/Figure5a.png", name(Graph) replace
graph export "$Figures/Figure5a.pdf", name(Graph) replace